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We present a model independent study of the chiral condensate evolution in a hadronic gas, in 
terms of temperature and baryon chemical potential. The meson-meson interactions are described 
within Chiral Perturbation Theory and the pion-nucleon interaction by means of Heavy Baryon Chi- 
ral Perturbation Theory, both at one loop, and nucleon-nucleon interactions can be safely neglected 
within our hadronic gas domain of validity. Together with the virial expansion, this provides a sys- 
tematic expansion at low temperatures and chemical potentials, which includes the physical quark 
masses. This can serve as a guideline for further studies on the lattice. We also obtain estimates 
of the critical line of temperature and chemical potential where the chiral condensate melts, which 
systematically lie somewhat higher than recent lattice calculations but are consistent with several 
hadronic models. We have also estimated uncertainties due to chiral parameters, heavier hadrons 
and higher orders through unitarized Chiral Perturbation Theory. 
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INTRODUCTION 

One of the most interesting questions about QCD is its 
phase diagram, and in particular, the transition from a 
hadron gas to a quark gluon plasma, and the restora- 
tion of chiral symmetry at finite baryon density. On 
the experimental side, there arc several experiments like 
CERES, KEK, STAR, SPS, LHC and RHIC that probe 
this transition at different values of temperature and 
baryon density. On the theoretical side, first princi- 
ple calculations on the quark gluon phase performed on 
the lattice have been traditionally hindered by the well 
known fermion determinant sign problem. The search for 
alternative approaches to overcome this problem, at least 
for some values of the chemical potential, have recently 
bolstered the activity in lattice QCD at finite density (see 
for example 0, 0, 13 Q and for a review and further 
references). However, lattice studies still present some 
difficulties in the infinite volume extrapolation and most 
importantly in that they do not use realistic values for 
the quark and hadron masses, particularly the pions and 
kaons, which are the lightest mesons and the most abun- 
dant at low temperatures. 

However, from the hadronic phase it is also possible to 
obtain model independent predictions by means of Chiral 
Perturbation Theory (ChPT) 0, 0, @ which is the low 
energy Effective Theory of QCD (see for introduc- 
tions and reviews). Let us recall that the spontaneous 
chiral symmetry breaking of QCD requires the existence 
of eight massless Goldstone Bosons that can be identified 
with the pions, kaons and the eta. Therefore they are the 
most relevant degrees of freedom at low energies. With 
these fields, the ChPT Lagrangian is built as the most 
general derivative expansion, over Aivf^ ~ f.2GeV (the 
symmetry breaking scale) , respecting the symmetry con- 
straints of QCD. Hence ChPT is the low energy effective 
theory of QCD. Actually, there is also an explicit symme- 



try breaking due to the small quark masses that give rise 
to a small mass for the pions, kaons and the eta, which 
are thus just pseudo Goldstone Bosons. For this reason 
ChPT is an expansion also in masses which are treated 
perturbatively. The leading term is fixed once we know 
the symmetry breaking scale, and the next orders contain 
a finite number of constants that absorb the infinities 
generated from loops, rendering the calculations finite 
order by order. Baryons can also be included in the La- 
grangian respecting chiral symmetry but their treatment 
is more involved due to their large masses. Within Heavy 
Baryon Chiral Perturbation Theory (HBChPT) [13 this 
problem is overcome for the meson-baryon interaction by 
an additional expansion over the baryon mass. Since we 
are only interested in the quark condensate in a hadronic 
gas, we will neglect the nucleon nucleon interactions. In 
particular, this leaves cold nuclear matter outside our 
domain of applicability. We therefore have a model in- 
dependent formalism derived from QCD that provides a 
systematic expansion at low energies and chemical po- 
tentials. This approach has proven very successful and 
works remarkably well within the meson sector, whereas 
for the meson-baryon sector, the convergence is some- 
what slower. 

In order to include the thermodynamic effects of the 
temperature and chemical potential, we will use another 
model independent approach, namely, the virial expan- 
sion [lfj, Il2l 1 1 3| . It is a simple and successful technique 
already applied to describe dilute gases made of inter- 
acting pions and other hadrons 12]. In contrast to 
lattice approaches, it is very straightforward to introduce 
the baryon chemical potential. For most thermal observ- 
ablcs it is enough to know the low energy scattering phase 
shifts of the particles within the gas, which could be taken 
from experiment, avoiding any model dependence. How- 
ever, since we are interested in the quark condensate, 
defined as a derivative of the pressure with respect to 
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the quark masses, one needs a model independent the- 
oretical description since it cannot be obtained directly 
from experiment. Thus, in this work we extend to finite 
baryon chemical potential, /i^, the model independent 
approach that combines the virial expansion and ChPT, 
a method already applied at = in 0,0>0- 

For the condensate, it is therefore particularly impor- 
tant to note that we are using the physical hadron masses 
in all our calculations. Let us also remark that the ex- 
pansion that relates hadron masses to quark masses, par- 
ticularly for pions, shows some of the best convergence 
properties in ChPT, much better than the energy expan- 
sion itself. We therefore hope that our results at low 
energy and chemical potential could serve as a guideline 
for a correct inclusion of mass effects in further studies, 
for instance, on the lattice, at least in the isospin limit. 

The plan of the work is as follows. In the next sec- 
tion we introduce in detail the virial formalism and make 
a rough estimation, using the free gas, of its applicabil- 
ity bounds in the temperature-chemical potential plane. 
With a brief introduction to ChPT and the hadron mass 
dependence on quark masses, we incorporate the inter- 
actions, which we do in two different scenarios. A pure 
SU(2) gas of pions and nucleons, where we include sys- 
tematically contributions up to a given order, and a more 
realistic gas including kaons, etas and heavier hadrons, 
where we neglect certain contributions due to the Boltz- 
mann suppression. We will show the size of the dif- 
ferent terms paying particular attention to the pion- 
nucleon interaction, and we will provide phenomenolog- 
ical parametrizations of the condensate melting temper- 
ature. In addition we will estimate by extrapolation the 
melting line in the temperature-chemical potential plane. 
Finally, we will calculate our uncertainties due to the im- 
perfect knowledge of chiral parameters and of the high 
energy behavior. The latter will be done by using uni- 
tarized chiral amplitudes. We will conclude with a sum- 
mary and discussion of our results. 



THE VIRIAL EXPANSION 

The thermodynamics of a system of hadrons is encoded 
in the grand canonical potential density z = €q — P, 
where we have explicitly eo, which is due to the exis- 
tence of a vacuum expectation value even at T = 0, and 
P stands for the pressure 0, H3| ■ In our case we are 
interested in a multicomponent relativistic gas made of 
pions, kaons, etas and nucleons. Later on we will intro- 
duce heavier hadrons. In addition, we will assume that 
only the strong interactions are relevant and that the 



expansion 



miHEi reads 



baryon density, defined as ns 



is conserved. For 



that reason, the pressure will depend on the temperature 
T and a baryon chemical potential fiB- Note that the 
non-strange quark chemical potential /i q — /is/3 is also 
frequently used in the literature. The relativistic virial 
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where £j = exp[/3(^ — Mi)], (3 = 1/T, Mi is the mass 
of the i species and /Zj = 0, ±/Ub for mesons, baryons 
(antibaryons), respectively. The coefficients 
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correspond simply to the virial expansion of the pressure 
for a free gas 



?P free = - J" ff / dpp* log 1 - THe-nVWf-rt 
i ^ Jo 

(3) 

where r]i = 1 for bosons and rji = — 1 for fermions. The 
interactions appear in the virial expansion through the S- 
matrix |l3l Il4j . For the meson-meson and pion-nucleon 
interactions relevant for this work, this can be recast in 
terms of the elastic scattering phase shifts, thus we can 
write: 

/3(Mi + M 3 ) r oc 

dEE 2 K 1 ((3E)A i:) (E), (4) 
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where K\ (x) is the first modified Bessel function of the 
second kind and: 

A«= £(2/ +l)(2J+l)Sf JS (E), (5) 

I,J,S 

Sy j s being the ij — > ij phase shifts (defined so that 
5 = at threshold) of the elastic scattering of a state 
ij with well defined isospin, total angular momentum 
and strangeness /, J, S. For the nucleon-nuclcon inter- 
action, using the 5-matrix representation could be more 
convenient than the phase shifts, but here we are only 
interested in a sufficiently diluted hadron gas, so that we 
can neglect NN, and the above formalism is enough for 
our purposes. Thus, in this work we will use the virial 
expansion together with ChPT phase shifts, extending 
previous works to include a non- vanishing baryon chem- 
ical potential in a hadron gas. We will study the validity 
of this approach in a section below. 

Let us now recall that the quark masses appear in 
the Lagrangian as m q qq, therefore the non-strange quark 
condensate is given by |14j : 

to)T» B = £ = <(%-#> - f£. (6) 

Note that we are working in the isospin limit, using a 
common mass m = (m u + rrid)^ both for the u and d 
quarks, and (0\qq\0) = (0\uu + dd\0). The (0\ss\0) con- 
densate, which is smaller, and whose thermal evolution is 
slower than the non strange one could also be stud- 
ied with similar methods, but this lies beyond our present 
scope. 
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ChPT 


IAM 






HBChPT 


HBChPT ± IAM 


"1 VP/ 


0.4 ± 0.3 


0.561 ± 0.008 


a i 




-2.60 ± 0.03 


-1.36 ± 0.02 


L r 2 (M p ) 


1.35 ±0.3 


1.21 ±0.001 


a 2 




1.40 ± 0.05 


0.438 ±0.015 


L:i 


-3.5 ± 1.1 


-2.79 ±0.02 


n-.i 




-1.00 ±0.06 


-0.70 ±0.04 


Ll{M p ) 


-0.3 ±0.5 


-0.36 ±0.02 


a 5 




3.30 ± 0.05 


1.29 ±0.04 


Ll(M p ) 


1.4 ± 0.5 


1.4 ± 0.02 


bi H 


- &2 


2.40 ±0.3 


3.06 ±0.3 


Ll{M p ) 


-0.2 ±0.3 


0.07 ±0.03 






-2.8 ±0.6 


-0.41 ±0.27 


L 7 


-0.4 ±0.2 


-0.44 ± 0.003 


b 6 




1.4 ±0.3 


-1.5 ±0.2 


Ll{M p ) 


0.9 ±0.3 


0.78 ±0.02 


bw 


-61s 


6.1 ±0.6 


7.4 ±0.5 


m(M p ) 


-3.4 ± 1.1 


-3.4 ± 1.1 


big 




-2.4 ±0.4 


-3.7 ±0.2 



TABLE I: One loop ChPT low energy constants used for our 
calculations. Those in the ChPT column are taken from Fiji] 
and , whereas those in the IAM column are taken from jl8j . 
HZ(M P ) is taken from ^ for both cases. 



TABLE II: Low energy constants for 0(p 3 ) HBChPT used in 
our calculations, for the non unitarized and unitarized cases. 
Taken from pcT ]. 



CHIRAL PERTURBATION THEORY 

Let us then briefly explain our use of Chiral Pertur- 
bation Theory (ChPT) and fix some notation, starting 
from the purely mesonic sector. Since the ChPT La- 
grangian is built as the most general derivative and mass 
expansion 0,0,0,0, the meson- meson interaction ampli- 
tudes are obtained as a series in even powers of momenta 
and masses, both denoted generically p 2 . The leading 
order starts at 0(p 2 ) and is universal, only depending 
on one scale F, which at leading order can be identi- 
fied with the pion decay constant /„-. At next to leading 
order (NLO) the amplitudes contain one-loop diagrams 
with 0(p 2 ) vertices plus the tree level contributions from 
the 0(p A ) Lagrangian. Indeed, all one-loop calculations 
can be renormalized in terms of a set of 0(p 4 ) parame- 
ters, Lk(/J.) and Hk(fJ-), fi being the renormalization scale, 
which can be determined from a few experiments, or in 
the case of , from the Resonance Saturation Hypoth- 
esis , and used for further predictions at low temper- 
atures. In Table I we list the values of these parameters 
that we will use for the amplitudes in this work. 

Within ChPT, partial waves tus are obtained as an 
expansion in even powers of momenta and masses. Drop- 
ping for simplicity the US indices, we find t(s) = 
his) ± £4(3) ± where t n (s) — 0(p n ). Since we need 
elastic amplitudes for the virial expansion, the phase shift 
5us(s) is simply the complex phase of its correspond- 
ing partial wave. In principle, the ChPT series is only 
valid at low energies compared with 47r/ w ~ 1 GeV, al- 
though in practice it is limited to momenta of the order 
of 200-300 MeV above threshold. For this reason we can- 
not apply our expansions for temperatures beyond that 
range. We will study the applicability of our approach 
with more detail in the next section. 

Baryons have a mass of the order of the chiral expan- 
sion scale 47r/,r, but they can also be included as degrees 
of freedom in a chiral effective Lagrangian if they are 
treated as heavy particles in a covariant framework called 



Heavy Baryon Chiral Perturbation Theory (HBChPT) 
[To). In this case, the series is organized in terms of order 
N = 1,2,3... that contain powers of p N / [F 21 M^ +1 -' 21 ), 
where I = 1, (N + l)/2. Once again, the divergences 
at each order are absorbed in a set of parameters, which 
in the notation of [2l| are called a's at 0{p 2 ) and b's at 
0(p 3 ). We provide in Table II the values we have used, 
but we want to remark that these parameters are strongly 
correlated and other sets can be found in [2^ |. 

The HBChPT framework is specially well suited for our 
approach, not only since ttN scattering has been calcu- 
lated to one loop, but also because it has been unitarized 
[2(3,|23| describing the A(1232) resonance correctly, which 
will be of use for the last section. In particular, we will 
use the third order calculation [2l|, |22j . Note that there 
is a fourth order calculation available |2^, that has also 
been unitarized |2oL |23| . but it introduces many more 
chiral parameters (up to 18 parameters in total) with 
very large uncertainties, and just an slight improvement 
over the third order results. Let us also remark that the 
second order coefficients, which are called cii in the for- 
malism of [^J that we follow here, can be translated to 
other coefficients often used in the literature, called Cj, 
which arc not dimcnsionlcss. The values listed in the 
HBChPT column in Table II correspond to: 

ci = -1.06 ±0.06 c 2 =3.4±0.5 
c 3 = -5.74 ±0.15 c 4 = 3.7±0.2, (7) 

in GeV -1 units. These are perfectly compatible with the 
recent values given in |26| 



ci 



C3 



-0.9; 
-4.7 



c 2 = 3.3±0.2 

1.2 „ _ o c+0.5 



-1.0 



C4 — 3.5_q 2 



(8) 



The corresponding values used for Unitarized HBChPT 
in the last section arc: 



ci 



C3 



-0.75 ±0.02 
-3.1 ±0.3 



C2 

c 4 = 



= 1.4 ±0.5 
1.5 ±0.2. 



(9) 
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As the last point concerning ChPT, and as we ad- 
vanced in the introduction, the chiral condensate virial 
calculation requires the knowledge of the quark mass de- 
pendence. Within ChPT, since we deal with hadrons, we 
actually recast Eq.© as follows 0| 

(«) = <o|^o>-£f^ (io) 

^— ' am oMh 

h 

where, for further convenience we have introduced the 
constant F, which is the pion decay constant in the chiral 
limit, and we have defined the coefficients 

BM 2 

c h = -F 2 ?-*(0\qq\0)-\ (12) 
am 

that encode the hadron mass dependence on the quark 
mass. The coefficient in SU(2) ChPT was calculated 
in 01 . Here we will use the SU(3) ChPT expressions 
at O(M^) which, for h = tt,K,t] can be read in [lfij |. 
Numerically they amount to 

c* = 0.9±g;t, c K = 0.5±° $, ^ = 0.418:?, (13) 

Using the O(Mj) calculation of the nucleon mass within 
HBChPT Hi Ej we find 

cat = 3.6±i;|, (14) 

where we have used the recent values [2f| for the 
HBChPT parameters in Eq.(JHJl. This corresponds to a 
nucleon sigma term: 

a nN = 40±2i MeV (15) 

where the uncertainties are estimated assuming uncorre- 
cted errors and provide a very conservative range. Had 
we used the constants in Eq.lQ, we would have found 
51^23 MeV. We have checked that the difference between 
these two central values will amount to roughly 1 MeV 
or less in our melting temperatures within the validity 
region of the approach, that we will study in the next 
section. 

Finally, let us analyze the hierarchy of the different 
terms. First of all, we see that the contributions to the 
virial expansion of the pressure are exponentially sup- 
pressed as exp(— 0Mh). Nevertheless, the exp(/3/z_e) fac- 
tor for the nucleons can overcome the previous thermal 
suppression if the chemical potential is of the same or- 
der as their mass. These exponentials are inherited by 
the derivatives. In addition, we observe that c^/M^ > 
cn/Mn » ck/Mk > c v /M v . For the above two rea- 
sons, pions and the interactions of the other species with 
pions, are the dominant contributions for studying the 
melting of the condensate. Nucleons become comparable 



only when their chemical potential is of the order of their 
mass. In addition, since nucleon-nucleon interactions are 
much stronger at low energies than those of pion-pion or 
pion-nucleon we will carefully exclude the region with too 
high nucleon density. 

As an illustration, and in order to compare with other 
results in the literature as well as to obtain an estimate on 
the reliability of the virial expansion, we will first study 
briefly the case of a free gas. 

THE FREE GAS AND THE DOMAIN OF 
VALIDITY OF OUR APPROACH 

The interest of this case is that we have a closed expres- 
sion for the pressure, already given in Eq. (J2J), which can 
be integrated numerically and compared with its second 
order virial expansion, namely, Eq. (|16fl with all £?™* = 0, 
i.e., 

pP^BPZh + BPtl (16) 

h 

First of all, in Figure 1 we have plotted the relative abun- 
dances of the most relevant species as a function of the 
temperature, which we parametrize in terms of their den- 
sity 

poo 2 

n h (n B ,T) = f^ dp P . (17) 

We see that the lightest particles, the pions, are the main 
component of the hadronic gas up to temperatures as 
high as 200 MeV. As previously remarked, we can also 
notice that nucleons form a small fraction of the gas un- 
less we reach high chemical potentials, of the order of 500 
MeV or higher. 

Next, in Figure 2 we show, in the (hb,T) plane the 
relative error between the exact calculation using Eq.J2J) 
and the second order virial expansion for the chiral con- 
densate, Ea. i|16|) . This we define as: 

exact (Qh) virial \ (18) 
ex act virial \ 

We see that the second order virial expansion provides 
a fairly good approximation e ~ 0.05 for moderate tem- 
peratures and chemical potentials (roughly T < 200 MeV 
and hb < 800 MeV) and slightly beyond we obtain just 
qualitatively correct results (e ~ 0.20), since the expan- 
sion deteriorates rather rapidly. In addition, and since 
our scope is just to describe a simple hadronic gas, the 
nuclear matter regime should be excluded from the valid- 
ity region. In fact, nuclear matter, even at low tempera- 
tures, is not a gas but a fluid 27], since, at low momenta, 
NN interactions are an order of magnitude larger than 
7T7T or TtN interactions. Just for illustration, we have plot- 
ted in Fig 2, as a continuous line, the points where the 
saturation density, po ~ 0.16 fm~ 3 , is reached. 



5 



0.5 



n, K, T) 

N ((ig = 40 MeV) 
N (u. B = 400 MeV) 
: 700 MeV) 




25 



50 75 100 125 150 
Temperature (MeV) 



175 200 



300 



250 "V ~ ■ 

> 200- 

I ; 

E 

3 150 - 

03 

o 

g 

£ 100- 

5Q 1 Validity region 



Po 
Po 
e = 
e = 



12 



0.05 
0.20 



■ \ \ 

. \ \ 
■ \ ' 
.\ 

\ \ 
'\ 

\ 

' "\ \ 

\ 



200 400 600 800 
Baryon chemical potential (MeV) 



1000 



FIG. 1: Density over T 3 of different species in a free hadronic 
gas. The pion population is by far the largest up to 200 
MeV, except at very high baryon chemical potentials, where 
nucleons dominate. 



Let us however remark that we are not interested in 
all the hadron thermodynamics, but just in a particular 
quantity, the quark condensate, that does not depend on 
the interactions themselves but on their derivative with 
respect to quark masses. Thus, the enhancement due to 
NN interactions with respect to irN or 7T7t is not as large 
for the quark condensate as it may be for other thermo- 
dynamic quantities that depend on the NN interactions. 
In particular, it has been shown [5^], using an equation 
of state for the nuclear matter regime based on a chiral 
Lagrangian, that the evolution of the quark condensate 
with nuclear density is well reproduced up to somewhat 
beyond po/2 just with a linear term proportional to the 
nuclear irN sigma term, but that higher orders in den- 
sity are relevant beyond (see Fig. 3 in |2jj). Since the 
nucleon-nucleon interaction lies outside the scope of this 
work, but the ttN interaction is included, we only con- 
sider densities below po/2. Note that with this choice we 
lie on the safe side. In addition, the uncertainties due 
to the irN interaction, that will be studied later in this 
work, are much bigger than the neglected effects of higher 
orders in density up to densities much larger than po/2. 

In summary, in Fig. 2, the dotted area stands for densi- 
ties above po/2 and the remaining white area in Fig. 2 will 
be referred as the "validity region" , in the understanding 
that this is just a very crude estimate. Fortunately, we 
will see that, once we introduce the interactions, most of 
the interesting phenomena, including our extrapolations 
for the condensate melting, occur within the bounds of 
this validity region, so that the virial expansion is reason- 
ably under control and definitely not diverging wildly. In 
particular, the region of relevance to study the freeze-out 



FIG. 2: Estimation of the region where the virial expansion 
may provide a good approximation to the quark condensate. 
Within the region inside the dashed line, and for a free gas 
of pions and nucleons, the relative error e of the second or- 
der virial expansion with respect to the analytic expression is 
less than 5 %. The virial approximation deteriorates rapidly 
beyond that "validity region". As commented in the text, we 
also require the nucleon density to be smaller than the sat- 
uration density po (the thick continuous line) and therefore 
we exclude the dotted area. The remaining white area is our 
approximate validity region 



or the phase transition in Relativistic Heavy Ion Colli- 
sions, hb = 40 - 50 MeV, T ~ 170 — 200 MeV H3 lies 
within our domain of validity. 

Finally, we show in Figure 3, the condensate melting 
line in the (/Lts,T) plane using either the closed form of 
the grand canonical potential density or the one extrapo- 
lated from our second order virial expansion. First of all, 
if we only consider n,K,r) and nucleons, we note that, 
starting from zero chemical potential, the second order 
virial melting line indeed follows closely the complete 
calculation with Eq.|J2J, although it deviates abruptly 
around T ~ 180 MeV and fi B = 800 MeV, that nev- 
ertheless lies beyond the "validity region" . 

In Figure 3 we also show the melting line that results 
if we add the contribution of the heavier hadrons in the 
free gas approximation by using 



A(gg) = - 



x ^ SAP dM h 

V dM h dm 

h 



(19) 



h 



dp ■ 



P 2 /th{p) 



■0(th{p)-Hh) 



Vh 



with e h {p) = vV + Mj. 

In such case, the extrapolated melting line falls slightly 
outside the estimated validity region. We will see in the 



Ch 


ChPT/HBChPT 


Nyffeler 


Tawfik-Toublan 


Cir 




2.0 


0.8 - 1.1 


CK 


U -°-0.7 


0.7 


1.6 - 2.1 


Cr, 




0.8 


2.4 - 3.2 


cn 


o fi +X.5 


4.2 


0.6-0.7 



300 r 



TABLE III: Coefficients Ch as calculated from estimates in 
other hadronic models |29l.l30| versus the analytic values ob- 
tained in ChPT or HBChPT. 



next sections that by adding interactions, an important 
part of the melting line moves within the validity region. 

Figure 3 is also relevant because it allows us to com- 
pare with previous results existing in the literature that 
use the free hadron gas, Eq.©, to study the condensate, 
including also heavier particles as above. In particular we 
can compare with the melting lines obtained in |2^| and 
pjof . We can observe that there is a quantitative differ- 
ence between those calculations and the one we present 
here. Such difference is mainly due to the dependence 
of the hadron masses on quark masses dMh/drh, which 
in those works was not obtained from ChPT, but sim- 
ply estimated as dMh/drh ~ N q (that is the number of 
light quarks, except for the pion that was chosen to be 
10) jig and 8M h /dMl ~ A/M h with A ~ 0.9 - 1.2 
|30j . Table 3 shows the comparison between the Ch co- 
efficients obtained from these estimations vs. the ChPT 
or HBChPT 0{p A ) calculations in Eas.lfT3)l and Ecis. ljTHl . 
which we consider very accurate, since the M^/(47r/ Tr ) 2 
ChPT expansion converges very well. Let us remark that 
the Ch coefficients in those works are generically larger 
than ours for h — n, K and rj. Hence, the effect of us- 
ing such estimates is to accelerate the thermal melting of 
the condensate. In particular, in Figure. 3, we show the 
results of |2^| , giving rise to a much lower critical temper- 
ature. This explains why the free gas results in [2(3, l3(l | 
yield melting temperatures much below those obtained 
in a free gas using ChPT calculations for dMh/drh, as in 
[lil llij. As a matter of fact, we have checked that our 
approach reproduces the results in |2^| just by changing 
the Ch coefficients. 

In what follows, and given the fact that the free and 
interacting contributions are separated in the virial ex- 
pansion, we show the results both expanding the free gas 
part to second order in the virial expansion, and also 
using the analytic expression. 



THE INTERACTING SU(2) GAS OF PIONS AND 
NUCLEONS 



In the literature it is frequently studied the SU(2) fla- 
vor case 0, EJ 0, E] , by considering the strange quark 
as heavy. This is simpler, since the pions are the only 
pseudo-Goldstone bosons and there is a clearer suppres- 
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FIG. 3: Condensate melting line in the (fiB,T) plane for 
a free pion and nucleon gas (continuous line). Note that 
it is very well approximated by the extrapolated line from 
the second order virial expansion (dashed) as long as fiB < 
800 MeV. The melting is faster when adding free heavier 
hadrons (dashed dotted line). Finally, we compare with the 
free gas results using crude estimates of dMh/dm q based on 
the number of valence quarks |2flj (dotted line). 



sion of heavier hadrons due to their heavy masses. On top 
of that SU(2) ChPT and HBChPT show a much better 
convergence than their SU(3) counterparts. In addition, 
we can use the one-loop calculations within HBChPT of 
the ttN scattering amplitudes. 

Thus, we will consider the second order virial expan- 
sion of a gas of pions and nucleons, where pions inter- 
act among themselves and with nucleons. The nucleon- 
nucleon contribution to the condensate melting is sup- 
pressed since c n /M w > cn/Mn and also because of the 
Boltzmann suppression of the nucleon population, which 
can only be compensated for baryon chemical potentials 
which lie outside the validity region of the virial expan- 
sion and are therefore beyond our scope. 

The hb = case within ChPT was first studied in 
|l4j | both with the virial expansion and by and effective 
field theory calculation of the grand canonical potential 
density. Later on the ChPT virial study was extended to 
a finite pion chemical potential in |15j . In both works, 
the only interacting particles were the pions, and all other 
hadrons were added in the free approximation. 

In this section, apart from the one loop ChPT 7T7t scat- 
tering amplitudes, we have included the 0(p 3 ) calcula- 
tion of nN within HBChPT H2- First we show in 
Fig. 4 the results at /is = 0, and we can notice how 
the free SU(2) gas melting temperature lies beyond the 
virial validity region, but the pion-pion interaction brings 
it down to ~ 230 MeV, within the naive validity region 
at hb = and in agreement with [Til Hjl fl6| . When we 
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FIG. 4: Chiral condensate thermal evolution at zero baryon 
chemical potential for the SU(2) pion and nucleon gas. The 
free gas extrapolated melting temperature falls entirely be- 
yond the validity region of the virial expansion, but 7T7t in- 
teractions bring it below 230 MeV. Although fis = 0, both 
the free nucleons and their interaction with pions accelerate 
sizably the condensate melting. Note that when interactions 
are introduced, the difference between calculating the free gas 
contributions exactly (resumed) or with the virial expansion 
is negligible, and the two curves fall on top of each other. 



was a difference of roughly 5 MeV observed at \ib = 
between the free gas melting temperature depending on 
whether we used the analytic form of the grand canoni- 
cal potential density or its second order virial expansion. 
The reason is that, in Fig. 3, after adding all interactions, 
the slope of the condensate is so steep that there is al- 
most no difference in the melting temperature, 1 MeV, 
between using the virial expansion for the free terms or 
not. 
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further introduce free nucleons we get an additional de- 
crease down to 224 MeV. Of course, within our approach, 
the melting temperatures are just extrapolations and lie 
near the edge of the validity region. What is not an ex- 
trapolation is the behavior at low T, where the different 
contributions are calculated consistently with the virial 
and chiral expansions. However, their effect is more diffi- 
cult to see in the figures, and we quote the extrapolated 
melting temperatures because it is easier to quantify the 
relative size of the different contributions. One should 
always keep in mind that these are just extrapolations. 

Surprisingly, the ttN interaction, which is the new con- 
tribution that we are adding, has a sizable effect even at 
fiB = 0, decreasing the melting temperature by another 
7 MeV, down to 217 MeV. 

Let us emphasize that each additional contribution de- 
creases further the melting temperature, but since the 
melting of the condensate accelerates near the melting 
point (the curve becomes steeper as it gets close to zero), 
even if the new contribution has a similar size as the pre- 
vious one, its effect on the melting temperature seems 
smaller. That is the reason why the relatively large ef- 
fect of ttN is even more surprising. A nice illustration 
of this effect is seen in Fig. 4 since the total result is 
practically the same no matter whether we use the virial 
expansion or the analytic (resumed) expression for the 
free terms. Indeed, both curves fall on top of each other 
in the figure. This is in contrast with Fig. 3, where there 



FIG. 5: Melting line in the ((j,b,T) plane of the chiral con- 
densate for an SU(2) gas of hadrons. We show the result for 
the free gas, but also adding the 7T7t ChPT SU(2) interac- 
tion to one loop, and the line resulting from adding the ttN 
interaction to third order in HBChPT. 

Let us now set /ij ^ 0. In Figure 5 we show in 
the (hb,T) plane the melting line of the free SU(2) 
gas of pions and nucleons, the gas with only interact- 
ing pions, and, finally, the effect of adding the ttN in- 
teraction. The effect of the latter becomes bigger as 
\xb increases from zero, and becomes maximum around 
[is = 600 — 700 MeV, where it produces a decrease in 
the melting temperature of about 40 MeV with respect 
to the gas without ttN interactions. Up to a chemical po- 
tential of (is = 40—50 MeV, the region of relevance |3jJ 
for Relativistic Heavy Ion Collisions, the decrease in the 
melting temperature amounts roughly to 10 MeV when 
we include the ttN interaction. 



In Fig. 6 we present a comparison of our pure 
results with lattice estimates in the literature 




Note that our melting temperature for a given \ib lies 
systematically above the corresponding one on the lat- 
tice. This could already be noticed in the /is = case 
[l4]| , both with the analytic ChPT calculation of the par- 
tition function and the virial expansion [bl. Il5| . 

To facilitate future comparison of our results with 
other calculations, we have performed a simple quadratic 
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tions that are numerically sizable. First, we consider the 
effect of free kaons and etas, which accelerate slightly 
the vanishing of the condensate, decreasing the extrap- 
olated melting temperature by ~ 4 MeV. Next, we have 
included the kaon and eta loops in the irir interaction. 
This amounts to calculate irir scattering in SU(3) ChPT 
instead of SU(2). Since the irir interaction is by far the 
dominant one, this one-loop effect decreases the melting 
temperature by ~ 2.5 MeV. When we include the one- 
loop interactions of ttK and 7r?y, we observe a further 
decrease of ~ 6.5 MeV. We have checked that, if we re- 
move the nucleons, these results are in agreement with 

To end with the \ib = , we add the heavier hadrons, 
where we have estimated that: 



dMh 
dm 



di 



(heavier hadrons) 



(21) 



FIG. 6: Comparison between our results and other estimates 
in the literature. The grey band corresponds to the fit in B, 
while the dashed line corresponds to the one shown in |32|. 
Note that these works obtain values for T c that are system- 
atically lower than our model independent approach (dotted 
line) . 



where a is an adimensional constant and d is the num- 
ber of valence u or d quarks in the hadron. For nucleons, 
following |3^|, and from our Eas. l|13|l and i|14fl . we esti- 
mate aw — 1-915 o- For the other heavier hadrons, we 
have estimated a ~ 0.5 — 2.5, which roughly covers the 
range of a obtained for the nucleons and light mesons. 



fit to our central results up to 500 MeV, giving the fol- 
lowing phenomenological expression for the dependence 
of the melting temperature with the baryon chemical po- 
tential: 

^M = l- 0.03155(2) f^V. (20) 
T c (0) V ' \T c {0)J V ' 

Where the uncertainty is only due to the statistical fit to 
the central value, and does not contain the uncertainties 
in our approach, mostly due to the chiral parameters, 
that will be treated in a separated section below. 

REALISTIC HADRON GAS 

In the previous section we considered a gas just made 
of pions and non-strange nucleons. However, in a real 
gas, we should consider all hadrons. This we will do 
by including them as free particles. The only exception 
will be the kaons and etas, that, in view of Fig. 1, are 
sufficiently abundant up to 200 MeV to deserve a sepa- 
rate treatment and include their interaction with a pion. 
All other interactions are severely suppressed by Boltz- 
mann and Ch/Mh factors. At very large \ib the heavier 
nucleons may not be Boltzmann suppressed, but that is 
beyond our estimated validity range. 

Thus, in Fig. 7, we compare the condensate thermal 
evolution at /j,b = of the SU(2) pion and nucleon 
gas including tttt and ttN interactions (thick continuous 
line) with the results obtained adding further contribu- 




Temperature (MeV) 



FIG. 7: Condensate thermal evolution at \ib — in a gas 
of hadrons. Starting from the SU(2) gas including nix and 
■nN interactions, we show the effect of adding kaons and etas, 
both as real and virtual states in nn loops, as well as the 
effect of the ttK and nw ChPT one loop interactions. The 
band includes the estimated uncertainties due to free heavier 
hadrons. 

Let us now turn to \xb 7^ . In Fig. 8 we show curves 
for the evolution of the chiral condensate for different 
constant values of the chemical potential. We can see 
how the increase of \xb accelerates the melting of the con- 
densate, and that the effect of the chemical potential on 
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the melting temperature is less than 10 MeV until /is is 
of order 250 MeV. In addition, we have plotted the curves 
as continuous lines within the "validity region" estimated 
in previous sections and as dotted lines outside. In this 
way we see that when /j,g is of the order of the nucleon 
mass, the aproach deteriorates rapidly, although it can 
still give reasonable results at very low temperatures. A 
different persrpective is provided in Fig. 9, where we now 
show the condensate evolution versus \i b for constant T. 
Note that, as we get close to T = the evolution with 
the chemical potential is almost a constant, the conden- 
sate in vaccuum, within the validity region. At exactly 
T = all fugacities vanish, except those of nucleons when 
Hb > Mn, but that lies beyond the validity region and 
cannot be studied with our approach. 
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FIG. 8: Condensate thermal evolution in a gas of hadrons 
for different constant values of /is. The dotted pieces of the 
curves lie beyond the "validity region" estimated in the text. 
Note the lines just correspond to central values. 



Finally, in Fig. 10, we present the melting lines in the 
(T, (is) plane, showing first the complete SU(2) pion and 
nucleon gas result, to which we have added all the effects 
of kaons and etas. As it could be expected, the biggest 
difference is observed at low \xb- In particular, in the 
interesting [lb — 40 — 50 MeV region [^H we find that 
there is a ~ 10 MeV additional decrease due to kaons and 
etas, on top of the one we already found in the previous 
section due to the ttN interaction. For higher energies 
the effect of K and ij are less evident not because they 
decrease, but because they become relatively less impor- 
tant with respect to the nucleon free terms and the due 
to the larger irN interaction. 

In the last section we will discuss the comparison of our 
results with some previous works in the literature, but 
first we will study the uncertainties within our approach. 



FIG. 9: Condensate evolution with the baryon chemical po- 
tential in a gas of hadrons for different constant values of T. 
The dotted pieces of the curves lie beyond the "validity re- 
gion" estimated in the text. Note the lines just correspond to 
central values. 
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FIG. 10: Our final result for the chiral condensate melting 
line in the (hb,T) plane. Starting from the SU(2) gas with 
7T7T interactions, we show the effect of irN interactions and of 
adding kaons and etas. The dark area covers the uncertainty 
due to heavier hadrons estimated as explained in the text. 



UNCERTAINTIES DUE TO THE CHIRAL 
PARAMETERS 



In previous sections we have shown results for the cen- 
tral values of the chiral parameters given in Table I and 
II. In this section we show the uncertainty due to the 
errors in the chiral parameters also listed in those Tables 
by adding in quadrature the uncertainties caused by each 
one of the parameters independently. Let us remark that 
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the dependence on these parameters is two-fold. On the 
one hand, they appear in the interactions through the 
phase shifts in Eq.JSJ, but, on the other hand, some of 
them also appear in the coefficients in Eq. <|12[) that 
parametrize the hadron mass dependence on the quark 
masses 16]. Indeed, the errors in Eas. (|T3|) and (CHEfl cor- 
respond to the uncertainties in the chiral parameters. 

Thus, in Fig. 11 we show again as a dark band our fi- 
nal extrapolated melting line in the (fJ.Bi T) plane, which 
includes the uncertainties due to heavier hadrons. The 
region between the dashed lines covers also the uncer- 
tainties due to the chiral parameters. The error is highly 
asymmetric since, as we have already commented, it be- 
comes more and more difficult to decrease the melting 
temperature by adding new effects, so that the inner 
dashed line is much closer to our "central" results. 
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FIG. 11: Our final result for the chiral condensate melting 
line in the {fiB , T) plane. The dark area covers the uncertain- 
ties due to the heavier hadrons, that we have included as free 
components of the gas. The area between dashed lines cor- 
responds to the uncertainty in the parameters of ChPT and 
HBChPT listed in Table I and II. 



that, as expected, and in contrast with lattice results, 
our melting temperatures always occur few MeV above 
the chemical freeze out temperature that has been esti- 
mated from p± spectra at RHIC as Tf = 165 ± 7 MeV, 
fig =41 ±5 MeV H3. 
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FIG. 12: Comparison between our results and others in the 
literature. Our melting temperatures are consistent within 
errors with ref. [3^| but systematically higher than ["113^1. 

Finally, to ease future comparisons with our results, 
we provide a simple quadratic fit to our results up to 
500 MeV, giving the following phenomenological expres- 
sion for the dependence of the melting temperature with 
the baryon chemical potential: 



Tc(0) 



1 - 0.02E 



h0.017 
-0.011 



fJ-B 

Tc(0) 



(22) 



In contrast with the fit given in the pure SU(2) case, 
Ea.Q2U|). the errors now are much larger since they take 
into account the uncertainty due to the chiral parameters 
and other heavier hadrons. 



We also remark that the uncertainties grow very fast at 
higher fis values in good agreement with the estimated 
range of validity of the chiral expansion. At sufficiently 
large values of /is our uncertainties are too large even 
for qualitative predictions. However, for small or mod- 
erate chemical potentials, say \ib << Mjv we get quite 
stable predictions, and the extrapolated melting line is 
determined up to a 10 to 15 MeV uncertainty. 

In Fig. 1121 we compare our melting line estimates with 
previous works. We show both our full error band (be- 
tween the dotted lines) and that due to heavier hadrons 
only (light grey band). Note that we are consistent 
within uncertainties with other hadronic models |34| . 
but systematically, our melting temperatures are some- 
what higher than estimates using lattice 0, Note 



UNITARIZED INTERACTIONS 

Up to this point we have been using a model indepen- 
dent approach. To that end we have used Chiral Per- 
turbation Theory, which is a low momenta expansion, to 
describe the interactions between hadrons. However, in 
the virial coefficients, Eq.Q the momentum integral ex- 
tends to infinity, and one may wonder if the extrapolation 
of ChPT amplitudes to high energies beyond their appli- 
cability range may be distorting the results. In order to 
check how important is the high energy contribution to 
those integrals we will thus extend the ChPT amplitudes 
by means of unitarization, and in particular the Inverse 
Amplitude Method (IAM), which provides a remarkably 
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good description of the meson-meson scattering data up 
to much higher energies. This will introduce a mild model 
dependency but will allow us to estimate how large is the 
effect of the high energy extrapolation in the virial inte- 
grals. 

Let us first briefly review unitarization and the IAM. 
For a detailed description of the method and its results we 
refer to 0, IH, US E3 • Note that, in view of Eqs.|@} and 
(0) , we need phase shifts of elastic amplitudes of definite 
isospin I and angular momentum J. Such amplitudes are 
called partial waves, tjj(s), although for brevity we will 
drop the /, J indices. Unitarity for elastic partial waves 
implies that 



Imt(s) = a(s)\t(s)\ 2 



a{s) = 2q CM /V~s, (23) 



where qcM is the center of mass momentum. Note that 
this implies that the modulus of the partial wave in the 
elastic regime is bounded by \t(s)\ < er -1 (s). 

The ChPT partial waves are obtained as an expansion 
in momenta and masses t(s) = t 2 (s) + ^(s) + where 
t 2 = 0(p 2 ), t 4 — 0(p A ), where p denotes a generic mo- 
menta or pseudo-Goldstone boson mass. Thus, unitarity 
is only satisfied order by order: 



lmi 2 (s)=0, Imk(s) =a(s)\t 2 (s)\ 2 , 



(24) 



This means that, since the truncated ChPT series grow 
like polynomials, unitarity can be badly violated at high 
energies. In addition, typical features of strong interac- 
tions, which are the resonances that saturate the above 
unitarity bound, are absent from the ChPT partial waves, 
since they correspond to poles in the second Riemann 
sheet that cannot be described with polynomials. In 
practice, pure NLO ChPT only provides a description of 
meson-meson scattering data up to energies of the order 
of 200 MeV above threshold. 

However, partial wave unitarity, Ea. (|23[l . fixes com- 
pletely the imaginary part of the inverse amplitude since 
it follows that Imi -1 = —a(s). Consequently, in the 
clastic regime a unitary partial wave can be written as 



t(s) 



1 



Reikis) -ia{s)' 



(25) 



and therefore we only need an approximation to 
Rei _1 (s). The Inverse Amplitude Method (IAM) sim- 
ply approximates Rei -1 (s) by the NLO ChPT result to 
find: 



t 



IAM 



(s) 



t 2 (s) 



t 2 — Re t 4 



iat\ 



l-U{s)/t 2 {sY 



(26) 



where in the last step we have used Ea. H24|) . In this 
way we satisfy unitarity while respecting the NLO chiral 
expansion, which is recovered at low energies. In ad- 
dition, the IAM can generate poles associated to reso- 
nances. The above derivation has been done for elas- 
tic amplitudes and yields amplitudes that indeed extend 



the ChPT calculations describing remarkably well the tttt 
and ttK scattering data up to roughly the first relevant 
inelastic threshold |35j . 

The formalism can be easily extended to the cou- 
pled channel case 0, |3(| H3- Let us illustrate the 
case when there are two coupled channel 1 and 2, (that 
could correspond, for instance, to 1 1) = \(tttt)u) and 
1 2) = \ (KK)[j)). We would then have four partial waves 
for the different choices of initial and final states tij , with 
i, j = 1, 2. Thus, the unitarity condition has now the fol- 
lowing matrix form: 



ImT — TT.T*, T 



t 2 l to 



a x 

(7 2 



(27) 

Following a similar argument as for the one channel case 
above, we find: 



T u = T 2 {T 2 -T A ) l T 2 , 



(28) 



This expression has been shown to provide a very good 
description of meson-meson scattering up to ^Js ~ 1.2 
GeV, much beyond the applicability range of standard 
ChPT. Note that this has been achieved with chiral pa- 
rameters Li compatible with those of ChPT and there- 
fore with a simultaneous description of the low energy 
and resonant regions [llil. In particular, the / (600), 
k(900), p(770), X*(892), / (980) and a (980) masses and 
widths are well described together with their associated 
poles in the second Riemann sheet [37) . 

Concerning the pion-nucleon interactions, the unita- 
rization could follow a similar line 23], but it can be 
improved [2(J taking into account that HBChPT has a 
double expansion: the low momenta chiral expansion, 
that counts powers of l/f^ and a heavy fermion expan- 
sion in powers of 1/Mm, where Mm is the nucleon mass. 
Factorizing explicitly the powers of Mm and / w , the ex- 
pansion now reads 



t = 



M«,(i.i) 
p 



Ml 



/-(1.2). 



N 



f 4 



N 



(29) 

f(ri,m) k em g dimensionless functions of uj/M^, where u> is 
the pion energy. Note that, contrary to the meson- meson 
case, odd orders of momenta do occur in the expansion. 

The partial wave unitarity condition with one fermion 
in the final state simply reads: Imt = qcAi\t\ 2 , so that 
now 



t(s) 



1 



Re t 1 (s) - iqcM' 



(30) 



Once again [23 we can use the HBChPT series at NLO 
to approximate Re£ -1 (s). However, the HBChPT series 
converges much worse than in the meson-meson sector, 
and for that reason it is convenient to reorder the heavy 
baryon expansion (see ref.[2(J), to get the following uni- 
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tarized amplitude: 



■IAM 



(31) 
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It has been shown that this expression greatly improves 
the plain HBChPT description, extending it at least to 
the nearest relevant inelastic threshold, and in particular 
up to y/s ~ 1400 MeV in the A(1232) resonance channel, 
which is the most relevant feature of irN scattering at 
low energies. Both the mass and width of this resonance 
are well described by the unitarized amplitude, together 
with its associated pole in the second Riemann sheet. 

We will now use Eas. (|28|) for meson- meson and (|31|) for 
ttN scattering, in order to extend ChPT at high energies 
and check the contribution of the high momenta region 
in the virial integrals. The chiral constants we used for 
the unitarized amplitudes are also listed in Tables I and 
II. Note that for meson-meson scattering they are com- 
patible with the values used in standard ChPT at NLO. 
The HBChPT parameters are not as well determined as 
those in the meson-meson sector, but it is important to 
note that the unitarization results are obtained with chi- 
ral parameters of natural size and quite compatible with 
different determinations from plain HBChPT. 

At this point we want to remark that since we are 
generating heavier resonances when using the unitarized 
amplitudes, we do not have to introduce them again as 
free particles in the gas, as we did in the previous sec- 
tion, to avoid double counting. In particular, we do not 
add now the / (600), k(900), p(770), if*(892), / (980), 
ao(980) and A(1232) resonances as free particles using 
Eci. (|20|l . since they are generated by the IAM. 

One advantage of unitarization is that some of the res- 
onances are now included with their actual width (that 
is, interacting with the mesons), which is an effect that 
we could not take into account by including them as 
free components of the gas. A priori, this could be 
relevant, since these resonances have strong decays, i.e. 
large widths, and considering them stable is one of the 
greatest over-simplifications in hadronic models. We are 
thus quantifying also that effect in our previous model 
independent parametrizations. Furthermore, we had a 
large uncertainty in their Ch or ah parameters, which 
is now largely reduced since they appear in the meson- 
meson or meson-nucleon virial terms, whose Ch coeffi- 
cients, Eqs. HI 3(1 and 1(14(1 . are much better determined. 

Thus, in Fig. 13, we compare the condensate thermal 
evolution at \ib =0 that we obtained in previous sections 
with the results using the unitarized interactions. The 
dark narrow band covers, for the IAM results, the un- 
certainty due to heavier hadrons for the non-unitarized 
results. As commented in the previous paragraph, the 
fact that unitarized interactions include the effect of the 



first heavier resonances, accelerates the melting. In ad- 
dition, since in the unitarized case what we call "heavier 
hadrons" are now the rest of hadrons, which are heavier 
and even less abundant, their associated uncertainty is 
much smaller. That is the reason why the dark band from 
the IAM is much narrower than the light band. Let us 
however, note that, once the uncertainties are included, 
the unitarized and non-unitarized results are very close, 
showing a remarkable stability of our results at /is = 0. 
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FIG. 13: Condensate thermal evolution at hb = with 
ChPT unitarized interactions. The resulting melting temper- 
ature is slightly lower but very close to the model independent 
non-unitarized case within uncertainties. 

Similar considerations hold for the uncertainty bands 
in Fig. 14. Again, the light band stands for the non unita- 
rized standard ChPT, and the much narrower dark band 
covers the uncertainty due to heavier hadrons of the IAM 
result. In this case, we see that the effect of unitarizing 
is larger, and the melting is systematically accelerated at 
moderate fi^ ^ 0. The bigger difference is now due to 
the fact that, as pointed out above, the HBChPT series 
converges much worse than the meson-meson one, and 
the effect of unitarization is more dramatic, particularly 
for the A(1232) that lies very close to threshold. Still 
these results give support to our statement that our melt- 
ing temperature estimates should be considered as upper 
bounds, although the differences between unitarized and 
non-unitarized are of about 10 MeV at most. 

Figure 15 shows the errors over this last band due to 
the uncertainties in the chiral parameters. The error is 
smaller than for the non-unitarized case. It should also 
be noted that both unitarized and non-unitarized cases 
overlap within their errors up to ~ 300 MeV and be- 
yond that are never further than two standard deviations 
away. 

Once again, the following quadratic fit provides a 
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FIG. 14: Condensate melting lines in the (/j,b,T) plane 
with unitarized ChPT interactions compared to the non- 
unitarized ChPT. The unitarized case melts faster than the 
non-unitarized, due to the dynamically generated resonances, 
that have their correct widths. 
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FIG. 15: The chiral condensate melting line in the (/is,T) 
plane for unitarized ChPT interactions. Note that the un- 
certainties due to the heavier hadrons (dark area) are much 
smaller than for the non-unitarized case in Fig. 11. The area 
between dashed lines corresponds to the uncertainty in the 
parameters of unitarized ChPT and HBChPT listed in Table 
I and II, and is also smaller than for the non-unitarized case. 



fairly good representation of our unitarized results up 
to 500 MeV 

Tc ^ B ) _ 1 _ n 032+ 0011 ( -J^-] (32) 
T c (0) 1 °- - ' 001 U(0)J ' ( j 

Where, once more, the errors take into account the uncer- 



tainty due to both the chiral parameters and the heavier 
hadrons 



SUMMARY AND DISCUSSION 

In this work we have presented a model independent 
description of a hadronic gas, based on the virial expan- 
sion and Chiral Perturbation Theory (ChPT) to NLO for 
the meson-meson interactions and third order for pion- 
nucleon within Heavy Baryon Chiral Perturbation The- 
ory (HBChPT). In particular, we have studied the melt- 
ing of the chiral condensate (qq) as a function of temper- 
ature, T, and baryon chemical potential, /ib- We have 
studied first the validity region, where the second order 
virial expansion is expected to provide a reasonably good 
description of the condensate evolution and nucleon den- 
sity is low enough so that we can neglect NN interac- 
tion. Since the virial expansion is a low density expansion 
and ChPT is a low energy effective theory, our approach 
is best applicable at low temperature and chemical po- 
tential, although it seems to provide still a fairly good 
description of the system up to temperatures of the or- 
der over T = 200 MeV at fj,B = and up to roughly 
fi B < 900 MeV at very low T (see Fig. 2). This allows 
us to quantify the size of different contributions to the 
evolution of the condensate in terms of its extrapolated 
melting temperature and chemical potential. Very likely 
this is a good estimate of the transition from "normal" 
hadronic matter, where chiral symmetry is spontaneously 
broken, to a phase where a chiral symmetry restoration 
may occur. 

We use the virial expansion because it allows us to 
include easily both the chemical potential and the in- 
teractions in the thermodynamics of the system. The 
chemical potential is very hard to implement on the lat- 
tice, to study rigorously QCD, and our results may give 
some guidance on this respect. The interactions are in- 
cluded from the zero temperature and density elastic 
phase shifts, which are fairly well known from experi- 
ment. As a matter of fact, for most thermodynamical 
quantities, one could simply take the phase shifts from 
data. However, the calculation of the condensate requires 
the knowledge of the dependence on the quark masses, 
which is not given by experiment. Hence, a theory is 
needed to describe those data, and since perturbative 
QCD itself does not describe the low energy regime, one 
must turn to ChPT, which is the low energy effective 
theory of QCD, for a model independent description. 

Our use of ChPT is very relevant since it provides 
model independent and very reliable calculations of the 
hadrons mass dependence on quark masses, which is a 
critical quantity to follow the condensate evolution. This 
dependence has been frequently estimated roughly in the 
literature, but the ChPT calculation provides an expan- 
sion in terms of meson masses which shows a good con- 
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vergence for pions, kaons, etas and nucleons, which are 
the dominant components of the hadron gas. Further- 
more, we have shown that several results obtained from 
free hadron gases in the literature can be reconciled with 
the interacting gas predictions if one uses these correct 
dependencies. 

Next, we have presented results for an SU(2) gas in- 
cluding only pions and nucleons. This is of interest be- 
cause its simplicity and because we can include system- 
atically all the effects up to a given order. This is in 
contrast with what we have called a realistic gas, where 
we have considered the interactions of kaons and etas 
with pions within the SU(3) ChPT formalism; although, 
formally, the KK, r\r\ and rjK interactions arc of the same 
order, we have neglected them due to: first, the thermal 
suppression of kaons and etas compared with the pions 
and, second, the weaker dependence of their mass with 
respect to non-strange quark mass. 

In this work we have shown explicitly the size of differ- 
ent contributions to the condensate melting for /_tg 7^ 0. 
In particular, we have evaluated the free kaon, eta and 
nucleon contributions, kaon and eta virtual effects in tttt 
interactions, as well as the irK, irrj interactions, these 
have also been calculated at fig = without nucleons 
to check with previous results in the literature. In addi- 
tion, we have included the ttN interaction, showing that 
it does have a sizable effect even at /ib = 0. Around 
[iB = 40 — 50 MeV the pion nucleon interaction itself de- 
creases the extrapolated melting temperature by 10 MeV. 
At higher baryon chemical potentials, the effect is larger. 

It is also relevant to remark that all contributions 
decrease the extrapolated melting temperatures. Intu- 
itively, adding more states increases the possibilities of 
generating disorder, i.e., entropy, and accelerates the 
melting of the condensate. The interactions seem to re- 
inforce this effect. Therefore, our extrapolated melting 
temperatures can be considered at least as upper bounds. 
Nevertheless, we have shown that it becomes harder and 
harder to decrease the melting temperatures with new 
contributions, so that we still consider that, within un- 
certainties, our results provide very good estimates of the 
melting temperatures. 

In this work we have also performed a detailed study of 
uncertainties, which are predominantly of two kinds: on 
the one hand we have approximated the effect of heavier 
hadrons by adding them to the gas as free components. 
The biggest uncertainty comes from the estimated de- 
pendence of their masses on the non-strange quark mass. 
On the other hand, we have large uncertainties from the 
chiral parameters of ChPT and HBChPT. This is the 
dominant source of error. All in all, our extrapolated 
melting temperatures come out with an error of the order 
of ±10 MeV, somewhat smaller at (Xb = and somewhat 
larger as /is increases. 

Let us remark that our results are consistent with other 
fj,B = model independent results that use ChPT and 



calculate the partition function either analytically |14| 
or also with the virial expansion 0, 0, 0] . We also 
find agreement within errors with hadronic models piij . 
but our melting temperatures come out systematically 
higher than estimates using lattice 0, . In addition, 
our melting curves lie, as expected, above the chemical 
freeze out temperature that has been estimated from p± 
spectra at RHIC jUJ. 

Finally, in order to have a more realistic description of 
the high momentum part of the virial integrals, we have 
used a unitarized version of ChPT that extends its appli- 
cability up to ~ 1.2 GeV, and of HBChPT that generates 
the A(1232) resonance. It also allows to describe cor- 
rectly the decay widths of the lightest resonances, which 
are the most abundant of the "heavier hadrons" that we 
had approximated as free components. We have found 
that our results are rather stable, within uncertainties, 
to this changes at lib = 0, but that these effects could 
produce a further decrease of the melting temperatures 
of ~ 10 MeV for higher chemical potentials. 

In conclusion, we have presented a model indepen- 
dent and systematic study of the chiral condensate evolu- 
tion with temperature and chemical potential, by means 
of the virial expansion and Chiral Perturbation Theory. 
The highlights of this approach with respect to other fun- 
damental approaches like QCD on the lattice are: the use 
of physical masses for the hadrons, the good control of 
the hadron mass dependence on quark masses, and the 
simple implementation of the baryon chemical potential. 
We hope that our results could serve as a guideline for 
future works in hadronic models and also on the lattice, 
in order to understand the phase diagram of QCD. 



NOTE ADDED 

After submitting to publication this work, we became 
aware of a recent lattice article 38], contested in 391, 
where they obtain at /_t B = 0, a T c = 192(7) (4) MeV, 
somewhat higher than all other lattice results we are 
aware of, and in fairly good agreement with the ChPT 
results [H, EH El] and here. 
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